New Consequences of Induced Transparency in a Double-A scheme: 
Destructive Interference In Four-wave Mixing 
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We investigate a four-state system interacting with long and short laser pulses in a weak probe 
beam approximation. We show that when all lasers are tuned to the exact unperturbed resonances, 
part of the four-wave mixing (FWM) field is strongly absorbed. The part which is not absorbed has 
the exact intensity required to destructively interfere with the excitation pathway involved in pro- 
ducing the FWM state. We show that with this three-photon destructive interference, the conversion 
efficiency can still be as high as 25%. Contrary to common belief, our calculation shows that this 
process, where an ideal one-photon electromagnetically induced transparency is established, is not 
most suitable for high efficiency conversion. With appropriate phase- matching and propagation dis- 
tance, and when the three-photon destructive interference does not occur, we show that the photon 
flux conversion efficiency is independent of probe intensity and can be close to 100%. In addition, 
we show clearly that the conversion efficiency is not determined by the maximum atomic coherence 
^ ' between two lower excited states, as commonly believed. It is the combination of phase-matching 

' and constructive interference involving the two terms arising in producing the mixing wave that is 

, the key element for the optimized FWM generation. Indeed, in this scheme no appreciable excited 

■ state is produced, so that the atomic coherence between states |0 > and |2 > is always very small. 

cn 
o 

CN ■ I. INTRODUCTION 

o . 

Efficient four-wave mixing (FWM) processes in the context of electromagnetically induced transparency (EIT) [1] 
Q^, involving a double- A scheme have been the subject of several recent studies [2-4]. In almost all of these studies, a 
large atomic coherence has been considered to be a key player in yielding a high conversion efficiency. Therefore, 
^ maximum atomic coherence was assumed. Furthermore, all studies available so far are based on either a steady-state 
2 approximation [2-4] or a full-numerical calculation[3] to obtain predictions on the system. A numerical calculation 
usually fails to yield as much insight into a problem as would be obtained from an analytical solution. The steady- 
• • . state treatment is analytical, but it oversimplifies the problem and, in some cases, it leads to incorrect conclusions. 
In particular, it is known that a steady-state treatment of a four-level double-A scheme can lead to an inconsistent 
solution in predicting quantum destructive interference effects that have profound relations to the wave mixing process. 
^ In this study, we describe an approximate analytical solution to a four-state double-A scheme that is different from 
5^ those studied previously [2-4]. Three features distinguish the present study from previous works: (1). A self-consistent 
fully time-dependent treatment leads to a three-photon destructive interference that extinguishes the laser excitation 
to the terminal state of a three-photon resonance, yet still provides nearly 25% frequency conversion efficiency. Such a 
three-photon destructive interference together with a high conversion efficiency from a quenched FWM state have not 
previously been reported in any double-A scheme in the context of EIT, and cannot be obtained from a simple steady- 
state treatment. (2). We show that resonance excitation (which leads to perfect one-photon EIT) is not the desired 
condition for optimum production of the mixing wave. We point out that sizable detunings from these resonances are 
required in order to avoid the three-photon destructive interference that limits the conversion efficiency, and (3) We 
show that, contrary to common beliefs, high conversion efficiency is not determined by having the maximum atomic 
coherence between the states [0 > and [1 >. Indeed, in this four-wave mixing scheme very little depletion of the 
ground state occurs and poi remains very small throughout the probe pulse. The correct criteria for achieving high 
conversion efficiency is a combination of the conventional phase-matching condition and a condition that ensures the 
constructive interference between two terms in the expression for the FWM field. In order to make both of the two 
terms of comparable size, a relatively large detuning is required to reduce near resonance absorption of one part of 
the four-wave mixing field. 



II. THEORETICAL MODEL 



Consider a four-level system that interacts with two transform-limited lasers {Eli{ujli) and El2{^L2), pulse length 
To, Fig.l). A short pulse probe laser (pulse length r << tq) is tuned near the [0 >— > [2 > resonance and fired at a 
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predetermined delay time. We assume that during the pulse length of the probe laser the two long pulsed lasers at ujli 
and ujl2 are sufficiently powerful to strongly saturate |2 >^ |1 > and |1 >^ |3 > transitions, therefore, |f2i2T| >> 1 
and |f2i3T| >> 1. As usual, flij = DijEL/{2h) and El are the one half Rabi frequency and the amplitude of the field 
for the respective transitions. Our objective is to seek a perturbative treatment for the response of the system to the 
short pulsed laser and to investigate the dynamics of the generated wave. 

We start with three equations of motion for the amplitudes of atomic wave function in the non-depleted ground 
state approximation. We thus require that |0o2| << \^2i\, so that most of the population remains in state |0 >. 
Taking Aq 1, we have 

— = ifli2A2 + ini3A3 + i{6i+i^)Ai, (la) 
dA2 



dt 
dA^ 



ii62 + h2/2)A2+in2o + i^2iAi, (lb) 
iiSa + 173/2)^3 + if23i^i + ^^^30• (Ic) 



at 

In Eqs.(l), Si is the detuning from the two-photon resonance between |0 > and |1 >, S2 is the dctiming of the short 
pulse probe laser from the |0 >^ |2 > resonance, 63 is the detuning from the three-photon resonance involving the 
transition |0 >— > |3 >. We will eventually assume that the lasers are tuned to the two-photon resonance |0 >— > |1 > 
and that the lifetime of this state is very long. Bear in mind that the generation of the four-wave mixing field 
will be shown to be very cfhcicnt, therefore, we must solve the above three equations of motion simultaneously 
with Maxwell's equations for both the probe and the generated fields. Taking unfocused beams and introducing 
Kij = 2TrojjiN\Dij\'^ / (he) where A'' is the concentration in cm~^, Maxwell's equations for these fields in slowly varying 
amplitude and phase approximation can be expressed as 

dfl20 1 dfl20 



dz c dt 



«Ko2^2, (2a) 



5^30 15^30 . , , 

Notice that in the present model, there exist two very different time scales, i.e. r << to, therefore, if the probe pulse 
occurs at the peak of the long pulse lasers, the amplitudes of the latter will remain nearly constant through the entire 
probe pulse. We therefore will be seeking a fully time dependent response of the system during the period when the 
probe pulse is present while treating both long pulse laser fields as time independent quantities. With this method in 
mind, Eqs.(l) and (2) can be solved analytically. Taking Fourier transform on the both sides of Eqs.(l) and (2), we 
obtain 

^21011 + {52+uj + z72/2)a2 = - W20, (3a) 
((5i -l-a; + i7i/2)ai + Oi2a2 + ^^130:3 = 0, (3b) 
9.3iai + {53+uj + i^3/2)a3 = -W30, (3c) 
dW2o .w . 

— 1—W20 = iKQ2a2, (3d) 

oz c 

dWso .u . 

— 1—W30 = iKosa^. (3e) 

oz c 

where ai, a2, a^, W20 and W30 are the Fourier transforms of Ai, A2, A3, ^20 and fiso, respectively. Eqs. (3a-3c) 
can be solved in terms of W20 and W30 with the result 

011 = ^— W20 ^— W30, (4a) 

012 = — W30 + W20, (4b) 

013 = — W20 + - W3Q, (4c) 

where, 

£>! = 5i -l-a; + 171/2, 

D2 = (^2 + 0.' + 172/2, 
D3 = 63+1^ + ^73/2, 
A = D1D2D3 - D3\^12\^ - D2\ni3\'' . 
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When these equations are used in Eqs (3d-3e), we obtain 

— 1-W20 = iKo2 7 W20-IK02 — T — 14^30, (5a) 

az c A A 

dWso .U . {\n^2?-DlD2) ,,, . ^31^12 ,,,, 

— i—Wio = IK03 T Wso-iKoa — T — W20- (5b) 

oz c A A 

For given W2o(0,a;) and with W3o(0,a;) = 0, Eqs (5) can be solved analytically, yielding 

Wso{z,u:) = ^l^^H^(^e^^^sin(Az), (6a) 
W2o{z,u;) = EH2^e^^^ (^^r^^-J^ sin(A^) + Acos(A^)) , (6b) 
where we have defined the new parameters 



A = ^(^^)' + aS3, (Ta) 
D = (7b) 

c A c A 



^^21^^13 r. ^^31^^12 
! ^ ,^3 - -K03 2 

If we use, for the probe laser at the entrance to the medium. 



'S'2 = -K02 ^ ,5'3 = -K03 ^ • (7d) 



Oo2(0,i) =no2(0,0)e-(*/^) , 

we find 

W^02(0,^?)/^!o2(0,0) = -^e-'''/^ 
where we have introduced the dimensionless variable r] = cot. We then get 



^30 (z, i) = 1^ r d,e-^-^e^-^^>e-^^iSM'-^, (8a) 

/A 1^20 (0,0) r , -nl . 
^^20 U, - = T=— / drie e 



iD{ri)z^-iri- 



.K2 - . . . . . 

I 1 sm(A^) + cosiA^;) 

2A 



In the following section, we will discuss the physical implications of Eq.(8) by examining some limiting cases where 
Eq.(8) can be evaluated analytically. In Section IV, we will evaluates Eqs. (8) numerically and compare the results 
with the analytical approximations. 

III. Discussions 

In this section we focus on the physical interpretation of Eqs. (8). We will consider some limiting cases where the 
inverse transform of the generated field can be carried out analytically. These limiting cases provide a great deal of 
insight into the wave propagation effects. 

Consider the limiting case where |rii2|, |f^i3| >> |(52|, |<^3|, 71, 72, 73- In fact, in what follows we shall always 
use 6\ = and assume that state 1 1 > is a second hyperfine level of the ground state, so that in an ultra-cold vapor 
7ir << 1. In this case the width of state |1 > is determined by the very slow rate of collisions between the atoms 
of cold low density vapor. Under these conditions, an accurate approximation to the parameters in Eqs. (7) can be 
obtained by assuming |Ai| << IA2I, IA3I, and |rii2|, |f^i3| >> IA2I, IA3I. Within this limit we can expand A by 
making use of the assumption that either |0i3|^, or |0i2|^ is much larger than either |£>iZ)2| or \D\Di\, therefore. 



^2+i^3 lf^2-^A .Qc^ n K12K32 

+ ^2*3 — ,7=; — To — 12- y^) 

«;i2|^13r + K32|f^l2r 
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If we replace sine and cosine in the inverse transform by complex exponentials, this approximation allows us to 
evaluate part of the integrals analytically for arbitrary pulse shape. When the detunings 62 and ^3 are very small 
compared with the half-Rabi frequencies, the other integrals are damped out in a small distance due to absorption. 

We first analyze the case where Si = S2 = S3 = 0, i.e. all lasers are tuned exactly on the unperturbed resonances. 
This is the limit where one-photon EIT is achieved. In this limit we find that after a large propagation distance, and 
at a point where the FWM has built up sufficiently 

^3o(-^, t) = |o''°fr^'io 12 ^20(0, t--^), (10a) 

^20{z,t) = ^03|»12P n^^^Qt-^)^ (10b) 



where 



1 1 _|_ K12K32 



Vgi C K12IO13P + AC32|r2i2P" 



This result indicates that in this limit ^^20(^5^) = {^2i/^3i)^3o{z,t).[5] Indeed, Eq. (4c) shows that in this limit 
a3 ~ for all ui, providing z is large enough to make | exp(i[Ko2|f^i3p + Ko3|f^i2p]/A])| << 1. It is required that 
I^^i2p >> I-D1-D2I in order for W02/VF03 = ^12/^13 to imply that 0:3 = 0. In other words, a destrucfAve interference 
has occured between the excitation pathways for the state. [6] We emphasize that this suppression of the excitation 
of the state |3 > cannot be obtained from the usual steady-state treatment of the atomic equations of motion. 
Indeed, the condition derived in the steady-state frame work for a three-photon destructive interference leads to an 
inconsistent prediction of the effect. We have evaluated the special case where K12T = K23T = 200cto~^, |r2i2T| = 5, 
|rii3T| = 20,7ir = 0.02, 72T = 73T — 2 using both Eq.(lO) and numerical integration of Eq.(8). As will be seen 
later the approximate analytical solution is in excellent agreement with the numerical evaluation and the ratio of 
fl2o{z, t)/V,3o{z, t) a.t z = 10 cm is equal to a constant (at all t where the two quantities are greater than 10~^ of their 
peak values) to an accuracy of seven significant figures. To this accuracy, the constant ratio is equal to Q,i2/fli3. Also, 
the peak in the two half-Rabi frequencies occurs at the point predicted by the group velocity Vgi . Note that if one 
allows |f2i3| = then Vgi reduces to the expression appropriate to the extremely slow wave propagation experiment 
[8], as should be the case. 

It should be pointed out that the destructive interference predicted above is equivalent to EIT for both the four-wave 
mixing photons and the probe laser photons. To see this clearly, consider the condition for Eq. (lb) to predict ^2 = 
at a point in space. This would be the case if Ai = — ri2o/^^2i at all times. With this amplitude for Ai the coupling 
terms dostriictively interfere. This is what happens when EIT occurs for the probe laser, for in this case there is no 
polarization to lowest order at cOp. Correspondingly, in order to have EIT at the four- wave mixing frequency we must 
have A3 = 0. Looking at Eq. (Ic), we see that this requires Ai = — O30/O31. In order for EIT to occur at both 
frequencies, the two values of Ai must be the same. This yields fl20 /^2i = ^30 /^3i • The same relation holds in Eqs. 
(10). 

We now consider the situation where |r2i2|T = |^^i3|t > 100, \Si\ = Oand7iT << 1. We also assume that \S3\T >> 1, 
but |(52/f^i2p, |<53/f^i2p, |<52/f^i3p, |i53/f^i3p < 1- These conditions indicate that both the Autler-Townes splittings 
and the detunings are large enough so that very little absorption occurs. In this limit three-photon destructive 
interference no longer occurs, and we find 

"3o(-,*) = .J'^ff'' ,2 (»2o(0,^-T^)-^2o(0,^-f)e^^-), (11a) 

'«02|iil3r + K03|i^l2r V Vgi Vg ) 

o ( ^\ K03|f^l2p 

^20(2, t) = 10 12' 10 12 (lit.) 



where we have introduced notations 



P 



i^l2|2<53 + |rii3P<52 



1 _ 1 ^ (l + |^^13/f^l2p)(«02|f^l3/f^l2p + K03) 



^9 C {S3^\^l3l^l2?S2) 



2 
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In the expression of Vg we have neglected a term that is of the same order of magnitude as 1/Vgi — 1/c, since 
such a term is much smaher than 1/Vg — 1/c in the above equation when the half-Rabi frequencies arc much larger 
than the detunings. Eq.(ll) indicates that there are two contributions to the growth of the FWM field fl^Q. The 
first contribution is due to the probe field that travels at the group velocity Vgi, whereas the second term consists 
of a probe field that travels at a second group velocity, Vg. If these two parts do not separate appreciably before a 
distance z such that \Pz\ = tt is reached, the two parts will interfere constructively. In the case where |r2i3/r2i2| — 1, 
this requires that c(ko2 + kos) << {^2 + ^3)^ and c/Vgi 1. The latter condition puts restrictions on kai2CT^ and 
'«32cr^, as compared with |ili2T| and ir^islr. 



The photon flux conversion efficiency under this condition can be found as 

(12) 



K03|^^12pK02|f^l3p ^_(^t-z/cf/T^_^iPz-Qz-{t-z/Vgf/T^^ 



Fp (K03|f^l2|2+K02|f^l3p)2 

where 

|ni2p73/2+|f]i3|272/2 



Q = P- 



11^1212^3 + |f^l3P^2 ■ 



In Eq.(12), Fm and Fp are the photon fluxes for the FWM and probe fields, respectively, and in deriving this 
expression we have, for the mathematical simplicity, assumed a Gaussian pulse shape for the probe field so that 
W2o{0,uj) = £^2ol|2)le-'^'^'. iProm Eq.(12), we notice that if we choose \Pz\ = tt and Ko2|f^i3p = Ko3|f^i2p, then 
we obtain F^ — Fp. This implies that the conversion efficiency is close to 100% whenever the difference between Vgi 
and Vg is small and the absorbtion factor \Qz\ << 1. Extensive numerical calculations have verified that conversion 
efficiencies of nearly 100% are indeed produced when the conditions stated before are well satisfied. 

In the case where the three-photon destructive interference does occur, i.e. in the case when perfect one-photon 
EIT is established with 5it = = SaT = 0, the second term in Eq.(12) disappears and the integral can be 
evaluated analytically. The result is a nearly 25% conversion efficiency in spite of the fact that the FWM field is 
strongly suppressed and there will be no further production of the field in the rest of medium. This indicates that for 
the scheme studied here, achieving one-photon EIT is detrimental to the high flux conversion efficiency. Therefore, 
one prefers to detune the two long pulsed lasers from the perspective resonances in order to avoid the destructive 
interference. Another important conclusion that can be immediately deduced from Eq.(9) is that it shows that 
maximum atomic coherence is not the optimum choice for maximum conversion efficiency as previously suggested [2] . 
Indeed, there are two conditions that must be met in order to achieve maximum conversion efficiency. The first one 
is the combination of k,q2, ^03, |i^i2|, and {flis] that maximize the amplitude of the expression in Eqs. (10-12). It has 
been shown [7] that this is precisely the conventional phase-matching condition for efficient wave mixing process. The 
second condition that must be met is that the propagation distance must be chosen properly in order to make the two 
terms in Eq.(12) interfere constructively. It is the combination of these two conditions that enables a highly efficient 
wave mixing and flux conversion process. 



III. NUMERICAL SIMULATIONS 



In this section, we investigate the model system numerically. Our main focus is to study the validity of Eq.(12) 
under specified conditions. This expression gives a complete description, under the conditions prescribed, on the 
flux conversion efficient with the field propagation effect included. It contains all the functional dependence of the 
FWM field on time, propagation distance, Rabi frequencies and detunings. We will first compare Eq.(12) with direct 
numerical evaluation of the inverse transform in order to show the validity of our approximation that leads to analytical 
results. After having established the validity of the analytical result, we then investigate the functional dependence 
of the conversion efficiency on parameters such as detunings and pumping Rabi frequencies, which may be useful in 
experimental verification of the theory. 

As mentioned before, the probe field propagation, generation of the FWM field, and the high-efficiency photon 
flux conversion described in the previous sections have all been vcrifled through extensive numerical calculations. In 
Figure 2a, we first show, for a typical set of parameters, the conversion efficiency predicted by the approximation 
(Eq.(12)) and by direct numerical integration of the inverse transform as a function of {t — z/c)/t. The parameters 
are so chosen that the three-photon destructive interference is defeated, therefore, both terms in Eq.(12) contribute 
to the overall conversion efficiency. Furthermore, we have chosen the parameters that maximizes the amplitude of 
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the expression given in Eq.(ll). The graph shows an excellent agreement between the approximate solution and the 
numerical evaluation of the inverse transform, and both methods predict nearly 100% conversion efficiency. In Figure 
2b, a different set of propagation parameters are chosen under the condition where the three-photon destructive 
interference is not present. Again, the plot shows an excellent agreement between the two methods. These results 
indicate that the analytical result Eq.(12) can, under the conditions given, correctly predict the propagation effect 
of the generated field. Therefore, we will use Eq.(12) to explore the functional dependence of the conversion efficient 
to various detunings and power densities. As has mentioned before, EIT process is established when both long pulse 
laser are tuned on resonance, i.e. = 5zt = 0. The resulting conversion efficiency can be analytically obtained. In 
Figure 3a, we show the dependence of the efficiency as function of S2T. As the detuning of the second laser from the 
resonance is increased, the efficiency increases, as expected. The similar effect due to the detuning d^r is also shown 
in Figure 3b. These figures indicate that for optimum operation, one should choose non-vanishing detunings which 
also satisfy the conditions forth set for the Rabi frequencies. Finally, we investigate the dependence of the conversion 
efficiency on laser powers. First, it is obvious that Eq.(9) is independent of the power of the probe laser. In Figure 
4, we plot the efficiency as function of |rii3/rii2p- As expected, the maximum conversion efficiency is achieved only 
at |r2i3/rii2p = 0.25, a necessary condition, for a given K02/K03 = 4, in order to achieve phase matching for the 
generated weave. 

V. Conclusion 

In conclusion, we have presented an approximate analytical solution to a four-level double-A scheme where laser 

induced transparency is expected. We predict that when the detunings from the three-photon resonance is small, 
and therefore the usual one-photon EIT process is established, the FWM field propagates without suffering any 
pulse distortion. Under these conditions. We also predict that at a sufficient propagation depth a complete three- 
photon destructive interference involving the FWM field and the three laser fields will reduce the polarization at 
the four-wave mixing frequency to zero. We have shown that even with such a robust three-photon destructive 
interference that strongly suppresses further generation of the FWM field beyond an onset propagation distance, the 
photon flux conversion efficiency can still be as high as 25%. Our result, however, indicates that one should detune 
both pump lasers from the perspective resonances, thereby avoid the usual one-photon EIT process. The use of a 
significant detuning will defeat a possible three-photon destructive interference (that limits the further production 
of the internally generated wave and therefore limits conversion efficiency to 25 %) therefore, allowing photon flux 
conversion efficiencies close to 100% with even a very weak probe field. In addition, our calculation shows that, 
contrary to common belief, a maximum atomic coherence between the two lower states is not the optimum condition 
for achieving maximum conversion efficiency. Indeed, with the current scheme \poi\ << 1. For high conversion 
efficiency, proper atomic parameters should be chosen according to Eqs.(ll) and (12). 
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Figure Captions 



Figure 1: Energy level diagram showing a four- level double- A scheme with relevant laser couplings. The theory and 
predictions presented in the text are expected to hold with minor modifications for other orderings of the energies of 
the excited states. 

Figure 2a: A plot of conversion efficiency comparing approximate solution and direct integration of the inverse 
transform as a function of {t — z/c)/t. Parameters used: |r^i2|T = 200, jriiajr = 100, kqsCt^ = 10, Ato2CT^ = 40, 
7iT — 72T = 73T = 0.1, SiT = 0, 62T — 20, S3T — 20, and z = 3.93 cm for constructive interference. The solid line is 
the approximate solution whereas the solid circles are full numerical solutions. Nearly 100% conversion efficiency is 
predicted by both methods. 

Figure 2b: Same plot as in Figure 2a except Sit = 0, S2T = 10, S^t = 10, and z = 1.96 cm for constructive interference. 
The solid line is the approximate solution whereas the solid circles are full numerical solutions. 

Figure 3a: Flux conversion efficiency vs. a dimensionless quantity 77 = (i — 2;/c)/t for a set of different detuning S2T. 
Open circle: S2T — 10, solid triangle: 62T = 20, open diamond: = 40, solid circle: 62T = 60. All other parameters 
are the same as in Figure 2a. 

Figure 3b: Flux conversion efficiency vs. a dimensionless quantity rj = (t — z/c)/t for a set of different detuning 63T. 
Open circle: (^sr = 10, solid triangle: S3T = 20, open diamond: S^t = 40, solid circle: 63T = 60. All other parameters 
are the same as in Figure 2a. 

Figure 4: Flux conversion efficiency vs. a dimensionless quantity rj = {t — z/c)/t for a set of different ratio of Rabi 
frequencies IO13/O12P. Cross: |ni3/ni2l^ = 1, open circle: IO13/O12P = 0.5, solid circle: |fii3/fii2p = 0.25, and 
open triangle: |fii3/fii2p = 0.1. All other parameters are the same as in Figure 2a. 
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